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Abstract 

The Two-Gluon Model of the Pomeron predicts strongly size-dependent high- 
energy hadron cross sections. Yet experimental cross sections for radially 
excited mesons appear surprisingly close in value. The strong coupling of these 
mesons in hadron collisions also predicted by the model permits a qualitative 
understanding of this puzzling behavior in terms of eigenmode propagation 
with a common eigen-cr. A detailed semiempirical coupled-channel model of 
the Pomeron is constructed to elucidate this and other features of high-energy 
hadron cross sections. 
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The Two-Gluon Model of the Pomeron (TGMP) has given a simple and rather successful 
picture of high-energy hadron-hadron scatterings Its most appealing feature is perhaps 
this: While accounting for the dependence on quark numbers emphasized so graphically by 
the additive quark model it goes beyond the latter by giving a natural explanation 

of the flavor dependence of hadron cross sections as a size effect arising from the color 
separation inside colorless hadrons This model prediction seems consistent with meson- 
nucleon (mN) cross sections for ground-state mesons |§. Theoretical support has come from 
nonperturbative QCD |§ and lattice calculations fl(| . 



The TGMP also predicts much larger mN cross sections for radially excited mesons of 
larger radii. In the small- meson limit 0, these cross sections are proportional to the meson 
ms radii, so that oyjv would be about four times larger than <T(j/^)n- 

Yet their experimental values appear to be similar: The older data on nuclear pho- 



toproduction, when analyzed in the Vector-Meson Dominance (VMD) model [jlT] gives a 
cross-section ratio of 0.77 (uncertainty 14) for the SLAC data [|12j, 0.77 (9) for the NA14 
data |TJJ, and 0.78 (11) for the E401 data |14j| . [When corrected with modern branching 



ratios [ 15] , these numbers should be changed to 0.79(15) and 0.86(10) for the first two data.] 



For muon production, one finds the ratio 0.81 (19) for the EMC data updated to 0.76 



'18) with modern BR's. For the E772 data on the proton production in nuclei [|T7|, analyzed 



in the Glauber multiple-scattering formalism [T8[ with uniform nucleon densities fitting ex- 
perimental nuclear radii [|l£j, we get 9.0 mb/7.8 mb = 1.15. A more careful analysis of this 
rough equality, already noted in |2(J , has recently been made in [^1[ . 



Yet another piece of the puzzle comes from the nuclear production amplitudes of excited 
mesons relative to that of their ground state in mN scatterings. Hufner and Kopeliovich 
(HK) [Q have recently deduced an experimental value of 0.46 for the production amplitude 
ratio ifj'/(J/ip) in p-nucleus (A) collisions |23|,|24|]. They have tried to explain this ratio in 
terms of a two-channel model of meson production using harmonic-oscillator (HO) wave 
functions and mN cross sections in the small-meson limit. We shall show below that this 
2-channel model is unstable, and that the generalization to many channels gives very poor 



results. 

The purpose of this paper is to show that this channel- coupling (CC) idea, when properly 
applied, seems to provide a key to a qualitative understanding of these puzzling features. 
The essential ingredients seem to be: (1) using the full TGMP without making the small- 
meson approximation, (2) using more realistic meson wave functions, and (3) using many 
meson channels. A semiempirical CC model of the Pomeron is constructed to elucidate the 
many interesting features of mN cross sections. 

Let us begin by pointing out that the coupling between different reaction channels of 
different hadron (h) excitations is unavoidable in hA scattering. What is not clear is whether 
the coupling is strong or weak in the TGMP. We shall see that it is strong for color-singlet 
mesons on nucleons, but weak for color-octet mesons. 

For definiteness, we consider the simplified situation where the target hadrons are the nu- 
cleons in a nucleus and where only the projectile meson undergoes radial excitations. Chan- 
nel coupling is in general a multiple-scattering phenomenon. Its physics becomes particularly 
transparent when described in the attenuation appoximation P5| , [Z2"| , where a multi-channel 
wave is attenuated by the matrix exp(— poLa) in channel space. Here po is the nucleon 
density, L = jR(l — 1/A) is the average attenuation length across a nucleus of radius R, 
and a is the matrix containing both single- and cross-channel absorption cross sections. It is 
usefully expressed as the dimensionless matrix M = <j/<jo, where o~ = cr n , so that M u = 1. 

The attenuation of the multi-channel hadron wave can then be understood in terms of 
the eigenvalues and eigenvectors of M. After a generally complicated transient, the wave will 
eventually decay into that eigenmode belonging to the smallest eigen-cr, hereafter called the 
"lowest" eigenmode. When produced in an eigenmode, the hadrons will propagate in that 
eigenmode over all distances, with the same eigen-a, and therefore the same attenuation, in 
all channels. However, only the lowest eigenmode is stable against perturbations. 

In the HK model ||22|| , the nucleons remain unexcited, but the propagating J/ip meson 



can be excited radially. In addition, the CC matrix M is calculated in the small-meson limit 
using HO wave functions. It is then a tridiagonal matrix in channel space. 



The HK model is based on the observation that in the 2-channel approximation, the 
theoretical ip f / ip ratio F of wave amplitudes in the lowest eigenmode agrees with the exper- 
imental value of 0.46 (6). This agreement is displayed in the first row of Table I, where Sj 
is the iih eigen-cx, and Vi(j) is the wave amplitude in the jth channel of its eigenvector. 

However, as the number n of channels increases beyond 2, the amplitude ratio F increases 
rapidly to above 1, thus destroying the agreement with experiment. At the same time, Si 
decreases towards zero, making the nucleus increasingly transparent to meson propagation. 
This result also contradicts the well-known experimental fact that the effective ipN cross 
section in nuclei is larger than its value in free space [^5[ . 

The table shows that in this HO r 2 model, stability in the ip' /ip amplitude ratio F requires 
the inclusion of at least 30 channels. The results numerically extrapolated to n = oo by 



using a "diagonal" rational approximant |26j are also given. 



The situation is more interesting if the small-meson approximation is not used. We start 
with the full hN amplitude of the TGMP first derived by Low and others [@,g§[/],|,|57|,[3B 
using either perturbative (P) or nonperturbative (NP or Cornwall) gluon propagators 



with an effective gluon mass m = m(q 2 = 0). Application of the optical theorem yields the 
total hN cross sections 

aij(h) = 8n h n N a 2 s J d 2 kD 2 (k)^(h; k)<S> N (k) , (1) 

where is the number of quarks in hadron i, and 

$ N (k) ~ 1 - f N (3k 2 ) , (2) 

^([qqhk) = S l] -f tJ (4k 2 ), (3) 

M»;*) = ^ + g/«(4fc 2 ). ( 4 ) 

M[(99)8y]i;*)^![^-/y(4fc 2 )]. (5) 
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Here is the diagonal or off-diagonal wave-function form factors. The proton and [qq]s 



expressions are from |[27|| , while that color factor for the hybrid is from pT| , |32| . The expression 
for hybrid mesons involves only the separation between g and (qq)s if the weak dependence 
on the qq separation is neglected. 

Each hadron factor in Eqs. (2-5) contains two terms: (1) a one-body term diagonal in 
the channel index caused by the exchange of both Pomeronic gluons with the same body: 
q, q, or the (qq)$ of the hybrid treated as a single body, and (2) a two-body term arising 
from the gluons interacting with both parts of the hadron, dependent on form factors, and 
responsible for channel coupling. For a colorless hadron, these terms interfere destructively 
because the scattering amplitude vanishes for a point hadron ||. This is why the cross 
sections are so sensitive to hadron sizes, and could be quite small even though each term is 
large. In color octets however, the two-body term is much weaker and adds constructively 
to the one-body term. The result is much larger cross sections and much weaker channel 
couplings. 

Ratios of these cross sections are theoretically simpler: They are independent of a s 
in the P treatment, and only weakly dependent on the QCD energy scale A.qcd m the 
NP treatment. We determine m by a best fit to the ratio irN/pN and KN/pN of the 
experimental Pomeron-exchange cross sections at the experimental hadron matter radii, both 
given in Table II. (Matter radii are deduced from charge radii.) The result is m = 0.08(10) 
GeV in the P treatment, and 0.27(6) GeV in the NP treatment used below, the latter being 



at the lower end of the Cornwall estimate of 0.5(3) GeV |29|. These masses are obtained 
with dipole form factors and Aq C d = 0.3 GeV, but the results for Gaussian form factors 
and other values of Aqcd are practically the same. Note that the product a s D(k) does not 
depend on a s in the NP treatment, and that the estimated error in m comes only from the 
errors in the input hadron radii, since the errors in the Pomeron cross sections are unknown. 

Meson wave functions are needed to calculate the CC matrix M. The effort is greatly 
reduced by using HO and the hydrogenic (Hy) wave functions. Neither turn out to be 
adequate, but they bracket more realistic wave functions, thereby allowing a simulation of 



the latter by the interpolation of their M matrices. As in Ref. ||22|| , we have used atotai 
instead of a a b s in M. 

As shown in Table I for J/ip mesons, the results for HO wave functions, though improved 
over the HO r 2 model, are not yet convergent with four channels. On the other hand, the 
results for Hy wave functions are beginning to converge with four channels, because the 
off-diagonal matrix elements are much weaker. This difference is also responsible for the 
result, already noticeable in Table I, that channel mixing in the lowest eigenmode becomes 
much weaker in the Hy model. (The channel components in the lowest eigenmode are highly 
coherent and have the same sign because all the off-diagonal matrix elements of M have the 
same negative sign.) 

The wave-amplitude ratio F for the first two channels turns out to be too large for the HO 
model, but too small for the Hy model, when compared to the experimental value of 0.46(6) 
quoted previously. To simulate more realistic wave functions, we take that linear combination 
of the two theoretical cross-section matrices which will reproduce the experimental value of 
F when the n = 2,3,4 results are numerically extrapolate to n = oo [J26|J . This gives the 
"75%(HO)/25%(Hy)" model shown in Table III. Obviously the fit would change with future 
changes in F, but the qualitative features of the model should be similar. 

The lowest eigen-a Si, which controls the steady-state propagation, should be lower 
than the smallest single-channel cross section ctq. It is about 0.64cr for J/ip mesons in the 
75/25 model. As channel-coupling decreases in going from the HO r 2 to the Hy model, the 
reduction of Ei below a also becomes less and less. 

The Pomeron-exchange parts <Jq of ground-state hN cross sections are themselves of 
interest. These can be obtained from the theoretical cross-section ratios by multiplication 
into the empirical pN Pomeron contribution of 21.70(s/GeV 2 ) a0808 mb fitted by Ref. ||. 
In this way, we not only reduce the considerable sensitivity of each cross section to model 
parameters, but also recover the experimental s dependence. Using dipole form factors 
fitting theoretical meson ms radii, we obtain the TGMP predictions shown in Table III. The 
errors shown come from the uncertainty of the gluon mass m. In addition, the cross sections 



are smaller by about 5% for Gaussian form factors. 

At y^s = 20 GeV, we find ctq — 5.2 mb for a J /if) meson ms radius of 0.044 fm 2 |]39 



This is about twice the value of 2.4 mb calculated in the gluon-fusion model of fl33|, and 



almost the same as the 6 mb obtained in ||25|| , or the 8 mb reported here, both from nuclear 
suppression data. If the meson ms radius is 0.055 fm 2 |40|| , ctq would have been 6.2 mb. 

Stable asymptotic propagation in nuclei on the other hand involves the lowest eigen-cr 
Ei. For J /if) mesons, this is ~ 0.64<To = 3.3(8) mb. The additional uncertainties from AF 
or the extrapolation to infinite matrix dimension probably does not exceed 0.3 mb each. 

A number of conclusions can be drawn: (1) Colorless mesons of different radial excitations 
and sizes experience the same total mN cross section in nuclear suppression if they are 
propagating close to an eigenmode, usually the lowest eigenmode. (2) The smallest eigen-cr 
Ei is smaller than cr = er n of the meson ground state. (3) For colorless J /if) mesons, Ei 
is smaller than the effective mN cross section deduced from nuclear suppression. (4) The 
single-channel cross section o"o is about twice the 2.4 mb given by the gluon-fusion model. 
Given the many uncertainties of our Pomeron model, it is not clear if the discrepancy is 
serious, and if so how it should be understood. Channel-coupling effects should also be 
present in the gluon-fusion model, but it is not known if they are strong or weak. (5) The 
nuclear propagation of hybrid mesons is quite similar to that for colorless mesons, but their 
cross sections are larger by the relative color factor |. On the other hand, color octets of 
different sizes tend to propagate independently in nuclei with roughly the same cross section. 
This is because the size-dependent, channel- coupling 2-body term is only -| of that in color 
singlets. 

Concerning conclusion (3), it has been suggested that the the larger cross sections of the 
absorption model could arise from the appearance of J / if) mesons in color octet forms, either 
together with a gluon in a hybrid ||32|| , or as bare color octets [53j]. A hybrid explanation of 



the increasing apparent nuclear absorption with increasing xf discussed in |54] can readily 
be given in our TGMP, with the ms hybrid radius increasing from about 0.03 fm 2 at small 
xf to about 0.25 (16) fm 2 at xf — 0.6. 



Conclusions (1) and (2) have interesting implications in the generalized VMD model 
of meson photoproduction from nuclei. In this model, the photon appears as a coherent 
admixture of all possible virtual vector mesons. Progress in this problem in the past has been 
hindered by the lack of information on both the channel-coupling matrix and the incident 
meson amplitudes |H |. Our CC model gives very specific predictions, although limited to 
only a small number of channels in the present calculation. The lowest eigenvectors of our 
largest n = 4 models, all with 75/25 mixing of wave functions, are shown in Table III for 
several meson families. The ground-state component vi(l) has been taken to be 1 for ease 
of comparison. Each eigenvector can be compared with the GVMD input amplitude vector 
v 7 = (1, /1//2, ■■■), where fi is the universal meson coupling constant to its source obtained 
from |15[ ) and for p', from . 



We see that the input amplitude vector for photoproduced mesons is rather close to the 
lowest eigenvector, especially for J/ip's. Hence mesons are produced close to this eigenmode, 
and propagate in nuclei with roughly the same reduced cross section Ex in different radial 
excitations. However, this reduction is not enough to explain why oyjmw is only 1.8-1.9 mb 
at y/s = 18-20 GeV when deduced from photoproduction data under traditional VMD |21[ 
unless a smaller gluon mass is used. 

The experimental ratio <j(ip') / 'cx(J /ip) for production in pA collisions is also independent 
of A ||23|| , meaning that the hadron production at a nucleon is also close to the lowest eigen- 
mode when interpreted in the coupled-channel model. This idea of "eigenmode production" , 
assumed both here and in [p^j , will required detailed justification. 

I would like to thank Dr. Cheuk-Yin Wong for many helpful discussions. 
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TABLES 

TABLE I. The lowest eigenmode of the channel-coupling matrix M for J/ip mesons in different 
models. 



Model 

n 


HO r 2 

Si/<7 


p* 


NP/HO 


p* 


NP/Hy 
Si/^o 


p* 


2 


0.61 


0.47 


0.66 


0.53 


0.92 


0.142 


3 


0.44 


0.68 


0.49 


0.72 


0.89 


0.153 


4 


0.35 


0.80 


0.40 


0.83 


0.88 


0.158 


15 


0.10 


1.10 










30 


0.05 


1.16 










oo 


0.00 


1.22 


0.02 


1.22 


0.86 


0.169 



*F = vi(2)/vi(l) 
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TABLE II. Pomeron-exchange contribution (in mb) to (ThN 
using dipole form factors. 



X hN s omm at ^ = 20 GeV 



Haclron 


(r 2 )* 


Ref 


XhN 




&hN 


D(k) 






NP 


P 


NP 



P 


0.67(2) 


[36 


7T 


0.44(1) 


[37 


K 


0.31(5) 


[38 


J/nj) 


0.044 


[39 




0.055 


[40 




0.181 


[39 


T 


0.013 


[39 




0.063 


[39 


p,u 


0.54 


[41 


K* 


0.37 


|41 


Ps 


0.54 


|41 




0.044 


[39 


T 8 


0.013 


[39 



21.70f 
13.63f 
11.82f 



35.21f 
22.12f 
19.18f 



3.2(7) 

3.8(8) 

8.6(9) 

1.1(3) 

4.2(9) 
14.90(12) 
12.6(4) 



3.7(5) 
4.3(6) 
8.8(7) 
1.6(3) 
4.7(6) 
14.91(9) 
12.6(3) 



5.2(12) 
6.2(13) 
14.0(15) 
1.9(6) 
6.8(14) 
24.2(2) 
20.5(7) 



28(8) 
29(8) 
29(8) 



45(12) 
48(13) 
47(13) 



In fm 2 . f Experimental results from ||. 
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TABLE III. The lowest eigenmode in the 75%(HO)/25%(Hy) model and the GVMD input 



vector v 7 


in different meson families. 


Meson ms radii used are those of Table II. 


Mesons 


Type 


Si/cr 


V! or v 7 


P 


n = 4 


0.69 


(1,0.69,0.50,0.35) 




n = oo 


0.46 


(1,0.90,...) 




7 




(1, 0.35, ...) 




n = 4 


0.71 


(1,0.40,0.19,0.09) 




n = oo 


n 6^ 






7 




(1,0.58,0.32,0.25,...) 


T 


n = 4 


0.74 


(1,0.32,0.12,0.05) 




n = oo 


0.70 


(1,0.35,...) 




7 




(1,0.64,0.58,0.40,0.45,0.29,...) 
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